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Abstract 

We compute eigenmodes of the lattice Dirac operator for quenched 
SU(3) gauge configurations on the 4-torus with topological charge ±1. 
We find a strong dependence of the zero modes on the boundary condi- 
tions which we impose for the Dirac operator. The lumps seen by the 
eigenmodes often change their position when changing the boundary 
conditions, while the local chirality of the lumps remains the same. 
Our results show that the zero mode of a charge ±1 configuration can 
couple to more than one object. We address the question whether 
these objects could be fractionally charged lumps. 



PACS: 11.15.Ha 

Key words: Lattice gauge theory, topology 



Introduction 

Understanding the nature of tiie QCD vacuum still is an open and chal- 
lenging problem. In particular, a possible connection between confinement 
and chiral symmetry breaking seems plausible. A hint that such a connec- 
tion might exist can be seen in the fact that at the QCD phase transition 
deconfinement and chiral symmetry restoration take place at the same tem- 
perature. 

A step in the direction of getting a grip on such a mechanism was the 
finding of so-called Kraan - van Baal (KvB) solutions (Ij. KvB solutions 
are self dual solutions of the classical Yang-Mills equations on a euclidean 
cylinder (S^ x IR'^). Thus they correspond to SU(3) gauge theory^ with 
temperature. A remarkable property of KvB solutions is that they consist 
of 3 sub- lumps and these three constituents together build up an object 
of topological charge 1. Furthermore these lumps were found to be BPS 
monopoles, a property that might lead to the possible connection to the 
confinement mechanism. Recently also KvB solutions with higher topolog- 
ical charge were found giving rise to 3\Q\ monopoles for the solutions 
with topological charge |Q|. 

Strong evidence that KvB solutions play a role also in the quantized the- 
ory was found by several studies on the lattice. For SU(2) analyses based 
on cooling with twisted boundary conditions ^3] and with regular boundary 
conditions ^ provided configurations that show the features of KvB solu- 
tions. The existence of constituents was established and the characteristic 
behavior of the Polyakov loop was seen. The zero modes of the Dirac oper- 
ator showed the expected change of location as a function of the boundary 
condition. For SU(3) evidence for KvB solutions comes from a study of 
eigenvectors of the Dirac operator El in thermalized (not cooled) config- 
urations. In j7j the zero mode of the Dirac operator in the background of a 
KvB soltion was computed and it was found that it strongly depends on the 
fermionic boundary condition used. The zero mode is localized on only one 
of the lumps but the particular constituent it chooses is determined by the 
boundary condition. This behavior of the zero mode was found for thermal- 
ized configurations in [HI and good agreement with the properties of KvB 
zero modes was established. It is furthermore known that the gauge lumps 
are predominantly (anti-) self dual [S] and a first step towards matching the 
results from cooling with fermionic methods is documented in 

The analysis of Dirac eigenmodes in thermalized configurations makes 

^Actually the KvB solutions are known for all SU(N), but here we only address the 
case of SU(3). 
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use of the fact that the low lying modes couple only to the long range 
structures of the gauge field and thus are an efficient filter for removing the 
UV fluctuations. Allowing for more general fermionic boundary conditions 
(in particular an arbitrary complex phase) makes this method even more 
powerful and led to establishing the KvB behavior of zero modes in finite 
temperature gauge configurations. 

It is natural to apply the generalized boundary condition technique also 
to gauge ensembles at zero temperature. A first step in this direction was 
made in [Hj where for a few configurations with topological charge ±1 it 
was shown that again the zero mode can be localized on different lumps 
when changing the boundary conditions. This hints at the existence of 
constituents also for charge ±1 configurations on the 4-torus T^. For this 
manifold no analytic solution is known (see ^01 for a discussions of the 
classical Yang-Mills fields on the 4-torus) and a lattice study has a more 
exploratory aspect in this case. However, for the case of two compact di- 
rections (T^ X m^) the recently found analytic solution [TJ again shows 
constituents and is already one step closer to the 4-torus analyzed here. 

In this article we present a large scale study of the Dirac eigenmodes for 
thermalized SU(3) configurations with topological charge ±1 on the 4-torus 
and analyze the dependence of the eigenmodes on the boundary conditions. 
We compare ensembles at two different volumes to analyze finite size effects. 
We show that in both volumes about 40% of the configurations have zero 
modes that change their position when switching from periodic to anti- 
periodic boundary conditions. This finding indicates that also on the torus 
a large portion of charge ±1 configurations gives rise to zero modes that 
couple to more than one object. 

An important question is whether these topological lumps have fractional 
charge and are truly constituents of a charge 1 object. The fact that the 
different lumps show up in the one and only zero mode we see supports this 
interpretation. Also the KvB solutions and a finite volume argument we 
discuss give hints in this direction. 

Description of observables 

For our calculations we use two ensembles of quenched gauge configurations 
generated with the Liischer-Weisz gauge action |12j on lattices with sizes 
12"^ and IG'*. For both lattices we use a coupling of /3 = 8.45 which gives rise 
to a lattice spacing of a = 0.094(1) as determined in (13) from the Sommer 
parameter. For the gauge fields the boundary conditions are periodic. 

We compute eigenvectors ip of the chirally improved lattice Dirac opera- 
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tor ^] which has good chiral properties and was shown [15] to reproduce the 
zero mode of a discretized instanton down to radii relatively small in lattice 
units. In our calculation we allow for a general phase at the temporal bound- 
ary of the lattice Dirac operator. The gauge fields have periodic boundary 
conditions and they do not single out a time direction so we choose the 
4-direction to be the time-direction where we impose the general boundary 
condition for the Dirac operator and its eigenmodes 

i;ix,U + l) = e^2-f ^^(^,1) • (1) 

The boundary phase parameter can assume values in the interval [0,1] 
with the two endpoints and 1 giving periodic boundary conditions. Anti- 
periodic temporal boundary conditions correspond to C = 0.5. 

For all gauge configurations generated we compute 50 (30 for 16^) eigen- 
values and eigenvectors with periodic boundary conditions {( = 0) using the 
Arnoldi method |16j . From the complete ensemble we then select those con- 
figurations which have exactly one zero-mode, i.e. have topological charge 
Q = ±1 according to the index theorem. For this subset of configurations 
we then solve the eigenvalue problem a second time now using anti-periodic 
boundary conditions (C = 0.5). The restriction to configurations with topo- 
logical charge Q = ±1 avoids possible problems with mixing of the zero 
modes and allows for a clean interpretation of the findings. Our statistics is 
291 configurations with charge Q = ±1 for 12^ and 199 configurations for 
16"^. For a small subsample we also analyzed the eigenvalue problems with 
increasing ( in small steps C = 0.0,0.1,0.2 ... 0.9. 

From the eigenmodes ip we can construct the scalar density p and the 
pseudoscalar density as 

p{x) = i'ix)a,ci^{x)a,c , (2) 

a,c 

P5{x) = '^ix)a,cil5)a,a''4^{x)a',c ■ (3) 

The star denotes complex conjugation. Both densities are gauge invariant 
observables. The scalar density p is always non-negative, while the pseu- 
doscalar density can have both signs, i.e. it is sensitive to the chirality of 
the lumps. Note that for the scalar density we have p{x) = 1; since the 
eigenvectors are normalized to 1. 

As we have discussed above, we find localized structures in the eigen- 
modes, i.e. localized lumps in p and p^. A first observable is the position 
^max q£ peak of the lump seen in the densities. As announced the 
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eigenmode can be located at different positions when comparing different 
boundary conditions. We will be particularly interested in the distance be- 
tween the maxima seen by the periodic and the anti-periodic zero mode 
which is given by 

I II max max n /a\ 

" II ■Aperiodic ~ •'^antip. 1 1 torus • \^) 

Here || .. || denotes the 4-dimensional euclidean distance. The subscript 
torus indicates that we always take the shortest distance possible on the 
4-torus. We also find that the eigenmodes come in different sizes ranging 
from tall and narrow (localized) to wide spread (delocalized). A convenient 
measure for the localization of an eigenmode is its inverse participation ratio 
/. It is defined as 

I = VY^pixf. (5) 

X 

Here V denotes the number of lattice points. It is easy to verify that for 
a maximally localized lump {p{x) = 5{x,xq)) one finds I = V, while for a 
maximally spread out lump {p{x) = 1/V) one obtains 1 = 1. Thus a large 
value of / corresponds to a localized, i.e. tall and narrow peak, while small 
values of / come from widely spread, delocalized modes. 



Properties of zero-modes with general boundary conditions 

Let us begin with describing some general features which we observed for 
the zero modes on the torus. We find that the zero modes show a single 
lump in both densities p and p^. For all values of C we find that the lumps 
are localized in space and time and on the average are rotationally invariant, 
i.e. no space-time direction is singled out by the shape of the lumps. 

When we change the boundary condition the number of zero modes 
remains invariant, i.e. we always find exactly 1 zero mode'^. When we now 
compare the zero modes with different boundary conditions but on the same 
configurations we often find drastic differences: The position of the lump 
in space time can be different and also the inverse participation ratio. In 
Fig. n we show a schematic plot of the typical behavior often observed when 
comparing the zero mode for periodic and anti-periodic boundary conditions. 

In the two plots of Fig. ^on the horizontal axes we show all of space time 
- the corresponding variable is x. On the vertical axes we plot some general 
density p. The zero mode is represented by the full curve and in this case the 
density p on the vertical axes is the scalar density defined in Eq. Q . The left 

■^We remark that for both lattice sizes we found for about 1% of the configurations an 
exception, i.e. a single value of ( where there was no or two zero modes. 
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Figure 1: Schematic sketch of the motion of the zero mode when switching 
the fermionic boundary condition from C = (l.h.s.) to Q = 0.5 (r.h.s.). 



hand side (l.h.s.) plot shows the case of periodic boundary conditions (C = 
0.0) while the r.h.s. plot is for anti-periodic boundary conditions (C = 0.5). 
For the different boundary conditions the zero mode is located at different 
positions in space time. When scanning through the boundary parameter 
( in small steps of 0.1 we found that the zero mode can be located at up 
to three different positions. To be more specific, for a subensemble of 10 
configurations on the 16^ lattice where we computed eigenvectors with small 
steps of C we found that 4 zero modes visited 3 different locations, 2 zero 
modes visited 2 locations and 4 zero modes remained at the same position 
for all values of (. 

For cooled configurations it has been established that lumps seen by the 
zero modes are associated with localized structures in the gauge fields. For 
thermalized configurations quantum fluctuations are present, but the zero 
modes and also the low lying modes couple only to the infrared structures 
in the gauge fleld JJj. In Fig. Q we indicate such localized lumps in the 
infrared part of the gauge field using dashed curves. Here the density on the 
vertical axis is e.g. the gluonic action density a quantity that is often used 
in the study of gauge lumps in cooled configurations. In the r.h.s. plot we 
also show the distance d (compare Eq. (jlj) between the lump seen by the 
zero mode with periodic boundary conditions and the anti-periodic lump. 

We stress, that when the zero mode visits different lumps when changing 
C, the lumps always have the same chirality. When inspecting the pseu- 
doscalar density we find that all lumps are seen with the same sign of p^. 
This is in agreement with the index theorem, which implies that if there is 
only one zero mode, its chirality is given by minus the topological charge. 
Since we do not touch the gauge configuration and only change (, the in- 
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dex theorem requires the chirahty of the zero mode to be the same for all 
fermionic boundary conditions. 

Let us discuss an important consistency check we performed. The bound- 
ary conditions for the gauge fields are periodic in all directions. In order 
to use the zero modes as an analyzing tool we single out one of the direc- 
tions by applying the boundary condition. However, the infrared structures 
of the gauge field should not depend on this choice. In order to check 
this, we applied the non-trivial fermionic boundary conditions not only in 
the 4-direction, but compared also with the zero-modes computed with the 
boundary condition applied to one of the other directions. We find that the 
same positions for the constituents are identified when applying the bound- 
ary conditions in different directions. 

Systematic comparison of periodic and anti-periodic zero modes 

Let us now come to presenting the results of a large scale comparison of 
anti-periodic and periodic zero modes. We begin with a discussion of the 
distribution of the distance d (compare Eq. Q)) between the peaks in the 
periodic and anti-periodic modes. The corresponding histograms for our 12^ 
ensemble are shown on the l.h.s. and for the 16^ ensemble on the r.h.s. of 
Fig. [21 The distance is measured in lattice units a = 0.094(1) fm and the 
histograms are normalized to 1. 

For both lattice sizes we find that the histogram has large contribu- 
tions for very small distances between d = and d = 2a. These contri- 
butions come from configurations where the periodic and the anti-periodic 
zero mode sit at essentially the same position, with the small discrepancies 
in the position partly due to quantum fluctuations. Clearly separated from 
these modes we find a well pronounced peak in the histogram containing 
configurations which show distances of at least d = 3a. For both volumes 
the number of configurations contributing to this peak is about half of the 
total statistics. 

An interesting question is whether one can understand the distribution 
of distances and the peak of the histograms. We experimented with a simple 
model where the positions of the two lumps are distributed independently 
on the lattice. In IR^ such a distribution would be proportional to but 
on the toroidal lattice we obtain a distribution which has a maximum at 
L/2 + 1 where L is the linear extension of the lattice. The corresponding 
curves are shown in the plots as connected dots. Note that these model 
curves are not normalized to 1. Instead their normalization was chosen 
such that they optimally describe (using a optimization) the data in 
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Figure 2: Distributions of the distance between the peaks in the zero modes 
seen with periodic and anti-periodic boundary conditions. The distance is 
given in lattice units. The l.h.s. plot is for lattice size 12^ (291 configura- 
tions), the r.h.s. plot for 16^ (199 configurations). The curve represents the 
numbers from a simple model to describe the data (see the text for details). 



the range d > 2a. Using this procedure we find that the area under the 
model curve is 0.425 for 12^ and 0.439 for 16^. For both volumes the model 
describes reasonably well the part of the histogram coming from separated 
periodic and anti-periodic zero modes. This implies that at least 40 percent 
of the charge ±1 configurations show two or more locations for their zero 
mode and the corresponding lumps seem to be distributed independently. 
Note that we only compared periodic and anti-periodic boundary conditions, 
and we thus do not count cases where the periodic and anti-periodic zero 
modes are located at the same position, but a different position is visited for 
some 7^ 0.0, 0.5. We thus expect that the fraction of configurations with 
"jumping" zero modes is actually larger than 40 percent. 

In Fig. we compare the distribution of the inverse participation ratio 
of the periodic zero modes to the distribution for the anti-periodic zero 
modes. The l.h.s. plot is for lattice size 12^, the r.h.s. plot for 16^. For both 
volumes we find that the distribution peaks at a relatively low value, with 
the distribution shifted to larger values of / for 16^. The remarkable feature 
of the plots is the good agreement of the distribution for the periodic zero 
modes with the distribution for the anti-periodic modes. This agreement 
demonstrates again that the size distribution of the lumps is governed by 
the infrared structures of the underlying gauge field and is not distorted by 
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Figure 3: Distributions of the inverse participation ratio for the zero modes. 
We compare the distribution for periodic (filled circles) and anti-periodic 
(open triangles) boundary conditions. 

the zero modes with different boundary conditions which we use as analyzing 
tool. Such a distortion could manifest itself by a different distribution of / 
for periodic and anti-periodic zero modes, which is clearly ruled out by our 
data. 

The nature of the topological lumps 

In the previous sections we have demonstrated that the zero mode can be 
located at different positions if one changes the fermionic boundary condi- 
tions. We have shown that for at least 40% of the configurations the max- 
imum of the scalar density visits two or more different space-time points. 
The zero modes are localized in all 4 directions and the chirality of the zero 
mode remains unchanged in accordance with the index theorem. The cru- 
cial question is now whether the zero modes hop between gauge lumps of 
fractional charge or simply move from one charge ±1 object to another. 

Let us first try to find a possible mechanism for the latter scenario, 
where the hopping comes about through some intricate interplay of inde- 
pendent topological objects with charge ±1. For example the appearance 
of three lumps with positive chirality seen by the zero mode could come 
from a configuration with three anti instantons and two instantons. Such a 
configuration would have net charge —1 and naively one would expect 5 zero 
modes. However, quantum fiuctuations could distort 4 of the zero modes, 
shift their eigenvalues up or down the imaginary axis and turn them into 
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near zero modes. The remaining single zero mode could then couple only to 
the three anti instantons and in this way show only positive chirality for all 
three lumps. 

This scenario cannot be entirely ruled out, but arguments were given in 
|18j that such behavior is unlikely to occur, at least so at finite temperature. 
Also our comparison of two different volumes makes this scenario unlikely. 
If one increases the volume, more topological objects with charge ±1 can fit 
into our box. This would imply that the scenario based on the interplay of 
charge ±1 objects could happen more often for the larger volume. Our data 
show however, that the rate remains close to 40% for both volumes although 
the larger lattice has a space-time volume more than three times as big as 
the smaller lattice. If on the other hand the hopping of the zero-modes 
is a genuine property of charge ±1 configurations, caused e.g. by fraction- 
ally charged constituents, then the percentage of charge ±1 configurations 
showing hops should be independent of the volume. This is what we observe. 

At this point we would like to comment on related work by Ilgenfritz, 
Martemyanov, Miiller-Preussker and Veselov I!) . In their article SU(2) 
gauge fields at varying temperature are subjected to a standard cooling 
technique. The authors concentrate on (approximate) solutions of the lat- 
tice field equations with topological charge ±1. Whereas in the confinement 
phase near the critical temperature a considerable fraction of configura- 
tions is seen to consist of separated constituents as static BPS monopoles, 
at lower temperatures topological lumps of integer topological charge, i.e. 
non-dissociated calorons show up. In terms of the Polyakov loop variable 
it turns out that these configurations have an intrinsic "±1 dipole" struc- 
ture connected to the nontrivial asymptotics of the holonomy. The form 
of the solutions is in one-to-one correspondence to the analytically known 
KvB solutions. In an analogous manner they have also studied the case 
of the 4-torus. They have found always non-dissociated "instantons" but 
with a similar dipole structure in all space-time directions for which the 
asymptotic holonomy is non-trivial. In this sense an interpretation of the 
topological lumps in terms of constituents is proven also in their paper. The 
fact that these lumps are not seen dissolved into constituents at T = is not 
straightforward to reconcile with our constituent interpretation of the zero- 
mode hopping for equilibrium fields. The lack of separate constituents in 
the cooled configurations may be an artifact of the standard cooling method 
which could drive the constituents into such a "bound state". Together 
with the authors of jl9j we are currently beginning a detailed analysis of 
this issue using APE blocking and improved cooling techniques combined 
with fermionic methods. 
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Summary and outlook 

In this article we show that for SU(3) gauge theory on the torus the zero 
mode of charge ±1 configurations can couple to diff'erent topological ob- 
jects. In particular we find that the zero mode often is localized at different 
positions in space-time when changing the boundary conditions, while the 
chirality of the lumps seen by the zero modes remains the same. For two dif- 
ferent volumes we compared the periodic and the anti-periodic zero modes 
for a large number of charge ±1 configurations. For both volumes we find 
that for about 40 % of the configurations the periodic and the anti-periodic 
zero modes are located at different positions. We compare the distribution 
of the distance between the different lumps to a simple model and find that 
the behavior of the constituents is compatible with independent placement 
on the torus. 

Although we cannot entirely rule out the possibility that the hopping is 
between integer charged objects, we lean towards the interpretation that it 
is true fractionally charged constituents we observe. This interpretation is 
supportet by analytic arguments at finite T showing that hopping between 
integer lumps is unlikely to occur. Also the absence of a volume dependence 
favors an explanation by fractionally charged constituents. 

Our results support a possible extension of the monopolc constituent 
picture known for manifolds with one or two compact dimensions to the 
4-torus. Such a generalization of the known classical solutions to the torus 
(with twisted boundary conditions) is certainly a challenging and interesting 
problem. On the numerical side also several interesting questions remain to 
be studied. In order to avoid problems with mixing of zero modes we have 
so far restricted ourselfes to configurations with exactly one zero mode. It 
would be interesting to analyze also the zero modes for higher topological 
sectors and check whether a similar behavior can be identified. Another 
direction for research would be an analysis of the near zero modes. They are 
expected to emerge when topological objects start to overlap. For the lowest 
modes one expects that the behavior of the lumps seen by the eigenvectors 
still resembles the properties of the exact zero modes. Showing constituent- 
like behavior for these near zero modes, which through the Banks-Casher 
relation are responsible for chiral symmetry breaking, would lead one step 
closer to a possible unification of the mechanisms for confinement and chiral 
symmetry breaking. 
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